Phase fluctuations, dissipation and superfluid stiffness in d-wave superconductors 
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o : 

_ We study the effect of dissipation on quantum and thermal phase fluctuations in 

£SJ ■ d-wave superconductors. Dissipation, arising from a nonzero low frequency optical 

_ ' conductivity which has been measured in experiments below T c , has two effects: (1) 

^\ a reduction of zero point phase fluctuations, and (2) a reduction of the temperature 

X^/^ ' a t which one crosses over to classical thermal fluctuations. For parameter values 

relevant to the cuprates, we show that the crossover temperature is still too large 
for classical phase fluctuations to play a significant role at low temperature. Quasi- 
particles are thus crucial in determining the linear temperature dependence of the 
i in-plane superfluid stiffness. Thermal phase fluctuations become important at higher 

■ temperatures and play a role near T c . 

o ■ 

O ; PACS numbers: 74.20.De, 74.72.-h, 74.25.Nf 

Oh; 

5a ■ I. INTRODUCTION 

c3 : 

_h \ There is considerable experimental evidence for a linear temperature dependence of the superfluid density at low 
Ch ■ temperatures in high-T c superconducting cuprates, i.e. 

Ps {T) = pM-aT (i) 

§ ' where a is a weakly doping-dependent constant However, there is still some controversy regarding the low- 

L ~ ~~ ', energy excitations responsible for this thermal suppression of p s . The simplest explanation is in terms of quasiparticle 

■ excitations near the e?-wave nodes |^-|(J. An alternative explanation is in terms of thermal fluctuations of the phase 
\ of the order parameter or other collective modes . 

The existence of well-defined quasiparticles in the superconducting state of cuprates is supported both by transport 
j 1 1 -|l3| and ARPES (w) experiments (even though there are some studies questioning their Fermi liquid description 
[ 15,|16|). However, the contributions of phase fluctuations to low temperature properties could still be important, 
especially in the underdoped regime where the superfluid density p s becomes vanishingly small as the Mott insulator 
is approached. 

The study of phase fluctuations raises several important issues: (1) The form of the phase-only action for layered 
d-wave superconductors (SC's) taking into account the long-range Coulomb interaction; (2) The crossover between 
quantum and classical regimes of phase fluctuations. These questions were studied in detail in ref. JTIj ] which, however, 
did not discuss the role of dissipation. In this paper, we focus on dissipative effects and how they affect the form of 
the phase-only action and the quantum-to-classical crossover. 

There are several reasons to believe that low energy dissipation is important in the high-T c cuprates. Theoretically, 
weak disorder within a self-consistent T-matrix calculation leads to a nonzero "universal" low frequency quasiparticle 
conductivity [ fl8| in <i-wave SC's. Experiments have also measured a nonzero low frequency conductivity p9| , ^6[ , much 
larger than the "universal" value. While there is no consensus on the origin of this large conductivity, one would 
definitely expect this dissipation to affect the phase fluctuations in the system, as first emphasized by Emery and 
Kivelson (EK) ||. Our formalism and results, however, differ from those of EK as discussed in detail in the paper. 
We summarize our main results below: 

(1) We derive the Gaussian effective action for phase fluctuations in the presence of dissipation using a functional in- 
tegral approach and integrating out fermionic degrees of freedom. While our effective action is derived microscopically 
by looking at fluctuations around a BCS mean field solution, we make contact with experiment by using parameter 
values relevant to the high T c SCs. We believe this phenomenological approach of using the derived form for the 
action, with coefficients taken from experiment, is valid for the SC state of the high T c materials, at least for T <C T c , 
when quasiparticles are well defined, fn addition, we also present in an Appendix, a hydrodynamic derivation for the 
phase mode based on a two-fluid model, which serves as a check on the microscopic derivation. 

(2) A dissipative quantum XY model is obtained by coarse-graining the Gaussian action to the scale of the coherence 
length and analyzed within a self-consistent harmonic approximation. 
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(3) We find that the magnitude of quantum fluctuations at T = is reduced by the presence of dissipation. 

(4) Dissipative phase fluctuations alone, in the absence of quasiparticle excitations, are shown to lead to a T 2 
reduction M of the superfluid stiffness. This behavior crosses over to a classical linear T reduction at a scale T c \, 
which decreases with increasing dissipation. 

(5) Choosing parameters appropriate to the high-T c cuprates, and overestimating the dissipation, we nevertheless 
find that the crossover scale T c \ is still fairly large. Thus one cannot attribute the low temperature linear reduction of 
p s (T) to classical phase fluctuations. This T-dependence must therefore arise entirely from quasiparticle excitations 
near the d-wave nodes. 



II. EFFECTIVE PHASE ACTION 

We find it convenient to express the superfluid density p s in terms of a stiffness D s — h 2 p s /m*, which in the 
London limit is related to the penetration depth A in a 3D bulk system through 1/A 2 = A-Ke 2 D s /h 2 c 2 . The in-plane 
superfluid stiffness in a layered system, with interlayer spacing d Cl is d c D s with dimension of energy. Henceforth, 
unless explicitly displayed, we set h — k B =1. 

We begin with the Gaussian phase action for a 3D isotropic superconductor (SC) 

Sg[0] = ^ J2 + D (^nW) 0(q, iWnM-q, -w n ). (2) 



For a derivation in the s-wave case see refs. |20|,[2lJ and for the d-wave case see ref. |17 . We show in Appendix A that 
the above action, and its generalization to layered systems, can be also derived from hydrodynamic considerations 
within a two-fluid model. In the above action (||), the compressibility x(q — * 0) ~ where is the Coulomb 

interaction and a is the lattice spacing. On continuing to real frequency D(u>) is the mean field stiffness, which is 
related to the mean-field complex conductivity a(u>) through D(u>) = (—iuja(uj)/e 2 ). In arriving at the above action, 
we have made the implicit assumption that <r(q, ui) ~ a(0,u>), and ignored the q-dependence of the conductivity |2^] 
for q < 7r/£ . 

We use the spectral representation for a, and find that 

D(iu n ) = D° s + i f° — ^ Rea reg (cv) (3) 

where we have used Recr(w) = irD®e 2 5(uj) +Hea reg (uj). For a frequency-independent Rea reg (u>) = <J DC , this simplifies 
to 

D(iu; n )=D s + ^f\oj n \. (4) 

Unless indicated otherwise, we will use this simplified form of the conductivity below, and use a = a DC d c /(e 2 /h) as 
a dimensionless measure of the dissipation. 

It is straightforward to generalize the above results to a layered system with an in-plane stiffness D® and a c-axis 
stiffness D° . Further, the Coulomb interaction in a system with layer spacing d c gets modified to p3f 



V(q) = 



2ire 2 d, 



sinh(gi|d c ) 



cosh(miG? c ) — cos{q±d c ) 



(5) 



where gii , q± are the in-plane and c-axis components of q respectively. 

To investigate the contribution of the phase fluctuations to the depletion of superfluid density it is necessary to 
go beyond the Gaussian approximation. The simplest model which allows for such an analysis is the quantum XY 
model, in which the phase field is defined on a coarse-grained lattice, with an in-plane lattice constant of £o and layer 
spacing d c . The coherence length enters as a short distance cutoff since the mean field assumption of a constant 
amplitude breaks down at shorter distances. 

Following exactly the same procedure of coarse-graining used in ref. |l7| (for the non-dissipative case) we now obtain 
the dissipative quantum XY action: 

^y[e] = ^^Y^^ + £k l | 7|| (Q)^(Q,c,„)| 2 + ^ J''* dr £ (1 - cos[0(R, r) - 0(R + a, r)]) 

D° d f £ \ 2 f 1/T 

- " ' ^" > ' x ' ' cos[0(R,t)-0(R + z,t)]). (6) 
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Here 7. (Q) = (4 — 2 cos Q x — 2cos<3 y ) with Q being the dimensionless momentum, and the scaled interaction ]24| ] 
V"(Q) = V(Q l] /£,o,Q±/dc)- While all momenta with \Q X \, \Q y \, |Q Z | < tt contribute in (||) above, the prime on the 
summation denotes a Matsubara frequency cutoff discussed below (see also ref. jL7] ) . 

In this derivation we have promoted the gradient terms in the Gaussian action arising from the superfluid stiffness 
to the cosine form, while the dissipative terms have still been retained at Gaussian level. A more sophisticated 
approach would probably end up with a t non-local kernel within the cosine term; we will however continue to work 
with the simplest action above. This action (^|) is well known in the literature as the resistively shunted Josephson 
junction (RSJJ) model and its phases and quantum phase transitions have been extensively studied [^5 26 1. Here we 
are interested in the effect of dissipation on quantum phase fluctuations and the classical crossover temperature, in 
the superconducting state. 

We now discuss the differences between our action (^) and that considered by Emery and Kivelson (EK) |S. EK 
included the effects of screening by replacing the V(Q) appearing in Eq. (||) with the screened interaction V S (Q) = 
V(Q)/e(uj). Here e(u>) = 1 + 47ri<7£(w)/a; is the dielectric function at Q = 0, and ctl(uj) is the longitudinal optical 
conductivity. Considering the isotropic 3D Coulomb interaction, 47re 2 /eocQ 2 , the dynamical term in the EK analysis 
reduces to the form w 2 eooe(w)Q 2 /47re 2 . This expression has been shown to be correct (27]]2lJ when the screening in 
the superconductor is produced by some external degrees of freedom with conductivity <jl(oj). An example is provided 
by a coupled system consisting of a superconductor interacting via Coulomb interactions with a normal metal. 

However, as discussed in ref. [^|, the EK effective action is not obtained for a single homogeneous SC. The 
longitudinal conductivity (7l of the SC does not explicitly appear in the expression of the density-density correlation 
function x- Instead dissipation appears through the transverse current-current correlation function, and affects the 
gradient term in the action (|^), so that D« — ► D, — iwcrj'(cj) / 'e 2 , where <tt is the transverse conductivity. 

If we would assume that physical (i.e., gauge-invariant) correlation functions appear as the coefficients in the phase 
action, then, using the equality of the physical longitudinal and transverse conductivities, it is easy to see that our 
action (§) is identical to the EK action ||, for the specific case of a SC with isotropic 3D Coulomb interactions. In 
this 3D case one could associate the a with either the gradient term, as we do, or with the time derivative term, as 
done by EK. The action used by EK is then formally the same as our action, and one could argue that dissipation 
should appear in the same way whether it is from an external bath (EK) or from internal degrees of freedom (our 
case). 

The above assumption of gauge invariant coefficients is, however, not valid for a single homogeneous SC where 
the screening arises from the (low energy) internal degrees of freedom. The coefficients in the phase action are then 
mean-field correlation functions; they cannot, in general, be gauge-invariant since the phase variable is yet to be 
integrated out. It is only upon integrating out the phase variable that one restores gauge invariance |17[|. 



III. QUANTUM AND CLASSICAL PHASE FLUCTUATIONS 



A. Variational Analysis 

We analyze the quantum XY action within the self consistent harmonic approximation (SCHA). We believe that 
this is adequate to calculate the effects of phase fluctuations at low temperatures, where longitudinal ( "spin- wave" ) 
fluctuations dominate and transverse (vortex) excitations are exponentially suppressed given their finite core energies. 
To examine the low temperature in-plane properties, we assume D Q ± = in (^) since it is very small in highly 
anisotropic systems with a large A^. For parameter values appropriate to i3z2212, we have numerically checked that 
setting = does not affect our in-plane results. 

The SCHA |2^,^Q is carried out by replacing the above action by a trial harmonic theory with the renormalized 
stiffness D.. chosen to minimize the free energy of the trial action. This leads to 

£>„ =^exp(-(^ 2 )/2) (7) 

where 59.. = (0 r _ T — # r + a , T ) with a — x,y and the expectation value evaluated in the renormalized harmonic theory 
is given by 

d 3 Q n " 7 (Q) 

(2^F n £ c ^ n e d c /V Q + (I?„ d c + % K|) 7|| (Q) " (8) 

As mentioned earlier, the dynamical phase distortions should have a frequency cutoff for the simple action we 
have considered. In our numerics, we use a cutoff fi~7j| n c corresponding to u) n < « / (£).. d c )(27re 2 /e 00 ^o)7 but we have 
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checked that the presence of a finite n c has only a minor quantitative effect on the results for {SO 2 ) in the presence of 
dissipation, and one may set n c — > oo to obtain qualitatively correct results. 

B. Analytical estimates of quantum and thermal fluctuations 

We first present estimates of the magnitude of quantum fluctuations and the thermal crossover scale making certain 
simplifying assumptions. The in-plane quantum fluctuations are seen to be dominated by relatively large from 

phase space considerations and the form of the integrand in (||). In this case, we may set V(Q ~ 1)/£q^c ~ 27re 2 /e 00 ^o- 
Restricting ourselves to low T, we ignore the Matsubara cutoffs and set n c — > oo. With these simplifications, we work 
in the limiting cases of small and large dissipation. We report further analytical results in Appendix B. In particular, 
we calculate the renormalization of the superfluid stiffness for an anisotropic 3D Coulomb interaction (instead of the 
Coulomb interaction in layered systems used in the paper) which permits us to analyze the case of arbitrary a. 

First recall the non-dissipative case Jl7[ | where the problem involves only two energy scales: the Coulomb energy 
(27re 2 /eoo£o) and the layer stiffness D..d c . The quantum zero point fluctuations of the phase are given by the dimen- 

sionless combination « / (2ire 2 / e^^o) / D d c , while the crossover to classical fluctuations takes place at a temperature 



T c \ ~ y/ (D^ d c )(27re 2 /e 00 ^o)- Taking into account the temperature dependence of the bare stiffness, a better estimate 
of the crossover temperature is T c \ ~ T c . 

For the dissipative case at T — we convert the Matsubara sum to an integral. For large a, ignoring £) d c in the 
integrand, and introducing a lower frequency cutoff, 27rD,, d c /a, it is easy to show that the the magnitude of quantum 

fluctuations may be estimated as (S9 2 ) — (=) In ^ yj '(27re 2 / 'eoo£o) / '(D y <i c ) . This is similar to the result obtained by 

Chakravarty et al J2(| for a RSJJ model with short range charging energies. Increasing dissipation thus leads to a 
decrease in quantum fluctuations as the system becomes more classical. 

To evaluate the temperature scale at which one crosses over to classical fluctuations in the presence of dissipation, 
we have to consider the temperature above which only the n = Matsubara frequency contributes, so that phase 
dynamics is unimportant. For large a, this may be estimated in a simple manner by setting (a/2ir)\uj n \ > D d c with 
n = 1, which ensures that fluctuations with n 3> 1 would contribute very little to the fluctuation integral in (0). This 
leads to T c \ > D..d c /a. A better estimate is obtained below, which gives T c \ « 3D^d c /W. It is clear that the classical 
limit emerges as the limit of infinite dissipation, a — > oo, for which T c \ — > 0. The crossover scale we obtain is similar 
in form to the estimate, T c \ ps T c /ef, given in ref. ||], but is much larger in magnitude. 

C. Low temperature behavior 

We next turn to the temperature dependence of the renormalized stiffness in the presence of dissipation, where 
we have set the bare stiffness to be independent of temperature. This is of course an unphysical assumption for a 
d-wave SC, but our aim is to explicitly check whether a linear T-dependence can be obtained within a model of purely 
dissipative phase fluctuations even when temperatures are smaller than the thermal crossover scale estimated above. 
The fluctuation (<56> 2 ) at low T can be evaluated analytically again by setting the cutoff n c — > oo. We can then cast 
the Matsubara sum in the form 

~ A(Q,T) ^ A(Q,T) A(Q,T) 

„_oo n2 + B (Q' T )H + C (Q' T ) ^ Q n 2 + B(Q,T)n + C(Q,T) C(Q,T)' U 

Rewriting the denominator of the first term in the form (n + ni)(n + nz), we separate out the terms using partial 
fractions and express the resulting sums in terms of digamma functions. As T — * 0, ni,2 — * oo which allows us to 
use the asymptotic expansion for the digamma function. The linear T term arising from the infinite sum is precisely 
canceled by the linear T term from the A(Q, T)/C(Q,T) term, leaving only a quadratic temperature dependence as 
was pointed out in ref. j§. We thus finally arrive at (58 2 )(T) = (56 2 )(0) + (V/3)(T/L> ||C ? C ) 2 at low T, from which 

^Kl-°(^—) 2 (10) 
D„(0) 6 [D^dcJ ■ { ) 

Thus, ignoring the effects of nodal quasiparticles, the asymptotic low temperature stiffness decreases as T 2 in the 
presence of dissipation. 
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At high temperature, above the thermal crossover scale, we recover the classical result (69*) (T) « 2T/(D || d c ). An 
improved estimate of the thermal crossover scale T c \ is obtained by matching the slope of this high temperature result 
for (S9^) with the low temperature result of Eq. (|ici|). This gives us T c \ = 3D ]f d c /a as stated earlier. 

At this stage, we turn to the recent results of Lemberger and co-workers |^9| , who use a circuit analogy and model a 
Josephson junction as an inductance (Lq) shunted by a resistance (i?) and capacitance (C). To make correspondence 
with this work, we note that the inductance Lq ~ (D dc)" 1 , the charging energy e 2 /2C ~ (e 2 /eb£o) and the resistance 
R ~ (l/<r)(/i/e 2 ). Upto numerical factors of order unity, our expressions for the magnitude of quantum fluctuations 
and the thermal crossover scale are then in agreement. The predicted quantum to thermal crossover has also been 
recently observed in experiments on conventional s-wave superconducting films pOfl, 



D. Numerical Results 



In order to obtain the various scales for the cuprates, we will choose parameters of the above action appropriate for 
the bilayer system Bi22\2 and evaluate the above estimates. We then present results of detailed numerical calculations 
which are shown to agree with these simple estimates. 

In the absence of detailed information on the bilayer couplings, we make the assumption that the two layers within 
a bilayer are strongly coupled and phase locked. Experimentally, the in-plane penetration depth of optimally doped 
Bi2212 is around 2100A and this translates into a bilayer stiffness ~ IhmeV. We use sw 10, and d c /a ss 4, with the 
in-plane coherence length £o/a sa 10. This leads to a Coulomb scale (e 2 /eoo£o) ~ 35meV. Using the above parameters, 
we find large quantum fluctuations in the non-dissipative case (ef = 0) with (<5# 2 ) ^ 1- The thermal crossover scale 
as estimated from the zero temperature bilayer stiffness, T c \ 3> T c ~ 100K . A more sensible estimate is obtained by 
considering the temperature dependence of the bare stiffness, and this leads to a crossover scale for T c \ ~ T c for W = 0. 

To study the effect of dissipation, we use conductivity data obtained from experiments performed in the supercon- 
ducting state. Consistent with our assumption of strongly coupled phases within a bilayer for Bi2212, the dissipation 
parameter a for this system will be taken to be the dimensionless bilayer conductivity. Recent measurements by Cor- 
son et al p6[ on Bi2212 films give a Drude conductivity with a large low frequency value corresponding to w 75 and 



a width of a few Terahertz. Similar large conductivities have been measured in the microwave regime |19 . We note 



that the "universal" quasiparticle conductivity predicted by Lee |lSj] for the bilayer conductivity (at T — 0,uj — » 0) 
corresponds to [g3| a m 24 for _Bz2212. The difference between this "universal" value, and the a 150 inferred from 
microwave data |19]l may be due to vertex corrections if^j . In our calculations, we use a constant dissipation with 
a ~ 150 (as an overestimate) over the entire frequency range of interest: oj n with \n\ < n c . This frequency range 
corresponds to oj < lOOmeV at T = 0. 

In the presence of dissipation, we find that quantum fluctuations are reduced to a very small value (<5(? 2 ) J$ 0.2. 
The thermal crossover scale is then T c \ = ZD„d c /a ~ 18K. This is consistent with our numerics, where we find that 
linear T behavior from thermal phase fluctuations only sets in above a temperature ~ 20K , for this magnitude of 
dissipation. While this is a low temperature scale, penetration depth measurements |j35| , fl9f observe a smooth linear 
T behavior down to much lower temperatures ~ 5K, which cannot be reconciled with this crossover scale. 

Turning to YBCO, and treating this system as weakly coupled single layers, far infra-red reflectance measurements 
pl[ appear to be consistent with a ~ 10-15 over a wide frequency range ~ 5-100mey. This smaller value of a 
compared to Bi2212 implies that dissipative effects are less important in YBCO. The "universal" conductivity |lj] 
for this case |p3fl c orresponds to a w 9, not inconsistent with the above data. However, low frequency microwave 
measurements [p2| on YBCO observe a strong frequency and temperature dependent quasiparticle conductivity, of 
the Drude form. We have checked that using a temperature dependent Drude conductivity in this very low frequency 
regime, in addition to a constant dissipation over the entire frequency range, does not significantly affect our results. 

In Fig. ^ we compare our analytical results and estimates obtained above, for i3z2212, with a numerical solution 
of the SCHA equations (Q) and (||) . We find that the purely quadratic dependence persists up to about 6K for the 
parameters discussed above, while linear T dependence only sets in at temperatures > 20K , consistent with the above 
estimate for T c \. Thus, it is impossible to ignore quasiparticles for understanding the smooth linear T behavior of the 
penetration depth which has been observed ]35| , [l9]| down to temperatures ~ 5K. 

We next include the linear T effect of quasiparticle excitations in the bare stiffness and ask how dissipative phase 
fluctuations renormalize this. The numerical result is plotted in Fig. || and shows that both the T = stiffness and 
its slope are renormalized by small amounts. This is completely consistent with our estimates for a small quantum 
renormalization in the presence of dissipation and a temperature scale of about 20K below which classical thermal 
effects are unimportant. 
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E. High temperature behavior 

Although quasiparticlcs dominate at low temperature, eventually phase fluctuations do become important at higher 
temperatures, in driving the transition to the non-superconducting state. The approximation used thus far (SCHA) by 
itself is clearly inadequate to address the problem of T c since it only includes longitudinal phase fluctuations. We thus 
proceed in two steps: first we calculate the temperature dependence of the superfluid stiffness within the SCHA, and 
then we supplement it with the Nelson-Kosterlitz condition for the universal jump in the stiffness at a 2D Kosterlitz- 
Thouless transition |3(| mediated by the unbinding of vortices. We do not take into account, for simplicity, the effect 
of layering in this calculation, but this could be easily done using well-known results for dimensional crossover. 

Our numerical results, obtained by solving the SCHA equations for parameter values relevant to i?i2212, are plotted 
in Fig. |[ We take the bare stiffness (the dashed line in Fig. ||) to decrease linearly with temperature, consistent with 
experiments, as a result of quasiparticle excitations. We extend this linear T dependence of the bare stiffness to 
higher temperatures assuming that the gap function is unaffected for T < T c . This is probably a good assumption in 
underdoped systems and may not be unreasonable for optimal doping. The renormalized stiffness calculated within 
SCHA is shown as the full line. This stiffness shows a jump at a temperature ~ 90K, but that is likely an artifact of 
the SCHA. The Kosterlitz-Thouless transition occurs at Tkt = ' K D s (T^ (T )d c /'& 1 where the fully renormalized stiffness 
D s is evaluated just below Tkt- We use the renormalized stiffness within the SCHA and the above condition, to 
obtain an approximate location of this transition. As can be seen from Fig. this gives us a reasonable estimate of 
T c sw 80K, when compared with experiments |35| , |l9fl on optimal Bi2212 which give T c sw 90K. 

IV. CONCLUSIONS 

In this paper we have derived and analyzed the effective action for phase fluctuations in the presence of dissipation 
arising from low-frequency absorption in a d-wave superconductor. We find that including dissipation reduces the 
magnitude of quantum phase fluctuations. The temperature at which one crosses over to thermal phase fluctuation 
is also reduced drastically. However, for parameter values relevant to the high-T c cuprates, we find that the thermal 
crossover scale is still large, so that quasiparticles dominate the asymptotic low temperature properties. In particular, 
they must be responsible for linear T-dependence of the low temperature penetration depth. 

Acknowledgements: We acknowledge C. Di Castro, M. Grilli, S. de Palo and T.V. Ramakrishnan for useful 
discussions and comments. We are particularly grateful to T. Lemberger for helpful discussions and detailed comments 
on the paper. The work of M.R. was supported in part by the D.S.T., Govt, of India, under the Swarnajayanti scheme. 
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APPENDIX A: TWO FLUID HYDRODYNAMICS AND COLLECTIVE MODES AT LOW 

TEMPERATURE 



In this Appendix we analyze the finite-temperature properties of the two-fluid model, to determine the collective 
modes of a superconductor in the presence of dissipative processes. We first consider a 3D Galilean invariant system 
and later generalize the result to a layered system, still maintaining Galilean invariance in the planes. Our goal is 
to determine the phase action that gives rise to this collective mode, which serves as another way to arrive at the 
Gaussian action derived microscopically in the text. 

The ordinary equations of a superfluid p7| are altered by the addition of a dissipative contribution. The longitudinal 
modes arising in the presence of long-range Coulomb forces obey the following set of linearized equations: 



u5p = q- j, 

ujSs + — q • (v s - v„) = 0, 
P 



. ep s 
mtu 



0Veg(w) 



iwv s — lqop H = 0, 

TO 

iq E = op. 



E 



(Al) 
(A2) 

(A3) 

(A4) 
(A5) 



Here the symbols have their usual meanings: w,q are the frequency and wavevector, the subscripts s,n refer to 
the superfluid or normal component, p, j are the particle density and particle-current density, v indicates a velocity, 
E is the internal longitudinal electric field, a reg (oj) is the regular part of the complex conductivity, and is the 
background dielectric constant. Further the thermodynamic variables, the pressure P, the entropy per particle s, and 
the chemical potential p, are related by the identity 



Sp 



-sST+-5P. 
P 



(A6) 



Note that dissipative processes due to thermal conductivity, which should appear on the right-hand side of Eq. (|A3|), 
can be neglected at low temperature and anyway, affect only second sound, which is decoupled from the density mode, 
as we shall see below. Moreover, the Lorentz force and the Joule heating are second order in the fluctuations, and do 
not affect the linearized equations which we are investigating. 

The equations for first and second sound, i.e., for density and entropy fluctuations respectively, are in principle 
coupled via the pressure variations: 



u> 2 5p - 



2c 2Ps s err 
U! OS = q 01 — 



2 , 



Aire p. 



5p + q 2 6P, 
4ir(jja reg (Lu) 



(A7) 



PsS x 

dp. 

PnP 



However, observing that we can rewrite 



op 



' dP 



ST, 



(A8) 



and using (dP/dT) p = as T — ► as a consequence of Nernst's theorem p^ ], we can conclude that, at low temperature, 



second sound does not mix with the longitudinal density modes. From equations (AS) and (A8), we then deduce the 
the dispersion for density fluctuations, given by 



47re 2 p s ,4iruj<7 re g(uj) 



i 2 2 

+ c„q 



(A9) 



where c p is a constant. From now on we neglect the term c 2 q 2 which is unimportant in the l ong wavelength limit. This 
is the plasmon dispersion, and coincides with result of the microscopic derivation in ref. fll7| ] . (To make connection 
with Eq. (29) of that reference, note that the real part of the right hand side of Eq. ( A9) above is ATrujlma(u!), where 
a(uj) is the total conductivity.) 
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Since we are interested in layered systems Eq. ( |A9| ) must be modified. For simplicity we consider the case of carriers 
confined to stacked (Galilean invariant) planes interacting via an the anisotropic Coulomb potential V(q) defined in 
Eq. (@). Referring to the components in-plane and across planes with a subscript || and _L respectively, we then 



assume that 



oo, mil = to, <7_l = and denoting a r « eg = an. The longitudinal electric field E has components, 



E\ 



-n\\9{% 



where the the electrostatic potential </>(q) for a density disturbance 5p(— q) is 
= oo and a± = 0, only the in-plane component Eu e nters in the li nea r resp onse 
response Eq. (A3) and determines the ballistic motion of the superffuid electrons, Eq. (A4). Rewriting (Al)-(Ra) and 



given by 0(q) 



and E± = -igj_0(q), 
= (V(q)/e)5p. As m± 



decoupling again the first and second sound as before, we obtain the density mode 



to e 2 



?fV(q). 



(A10) 



Eq. (AlO) allows us to deduce the correct expression for the phase-only Lagrangian in the layered case. Indeed, 
according to the previous equation, in the presence of dissipation, the superfluid stiffness must be transformed as 
Ps/rn — > [Ps/m — ioja\\(u})/e 2 ) . We thus obtain, at the Gaussian level, the Lagrangian density 



£(q,w) = o 



Ps (w) 



l^(q^)l 2 , 



(All) 



which is the same as the Gaussian action used in the paper. 



APPENDIX B: ANALYTICAL RESULTS FOR PHASE FLUCTUATIONS WITH ARBITRARY a 



It is useful to calculate the the fluctuation (<5# 2 ) for arbitrary dissipation er, even approximately. To make progress 
we work with the anisotropic Coulomb interaction 



4ne 2 



e u q 2 + s ± q 2 



(Bl) 



which is differe nt f rom the Coulomb potential (^|) for layered systems used in the paper. In particular, one cannot take 
a 2D limit of (Bl) unlike in the layered case. However, it permits us a simple evaluation of the integrals appearing 
for (68 2 ) in (|). 

On using the appropriate scaled Coulomb interaction and simplifying 7||(Q) — Q 2 , the renormalized Gaussian 
action (with D ± = 0) takes the form 



8T 



Aire 2 



Q 2 



D \\ d c + Trl^n 

Z7T 



Q 2 MQ,n) 



(B2) 



with rj = [e ± ^ 2 ) / {£ ^d 2 ) . Setting uj c = Aire 2 /d c , we then obtain the fluctuation 
(66 2 ) = 2Toj 



* dQ ± d 2 Q^ v 



£,i (1 + vCq)^u + [ D^d c + —\cu r , 



(B3) 



with c Q = Q 2 JQ 2 . 

The integrand in (|B3| ) depends on Q only through (q . The Q integral can therefore be transformed to an integral 
over the variable C, with density N(Q = 1/3VC for C < 1 and N(Q = l/3( 2 for C > 1, and 



) = / dCiV(C)*((l + V(h\) = / d(N(0 = $(e), 



where 



$(e) = 2Tu; c ^ 



zwz + ujr D„ d, 



2vr' 



-, -1 



(B4) 



(B5) 
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and e is a suitable average value of (1 + ryC) £ N ■ We have written (B4) such that the effect of the Q integral appears 
as a renormalization of the bare in-plane dielectric constant en to a larger value e. Of course, the value of e is 
temperature dependent. However, the leading temperature dependence of (50 2 ) is, in most cases, e-independent (with 
the noticeable exception of the dissipation-less case a = 0). Therefore, to proceed further analytically, we take e to 
be a constant henceforth. 

The quantum corrections can be calculated analytically, expressing (59 2 ) by means of the spectral representation 



for the Matsubara phase propagator deduced from Eq. (B5), so that 



( < >(T) = -L M Z ^yf^N^) 2 COt ^) ' (B6) 

where ui 2 — {D^d c ) (uj c /e) = Aire 2 p s /em* . 

The above expression allows us to explicitly calculate the value of (59 2 )(T = 0). Introducing the dimensionless 
parameter s = (er/47r) (lo c /£lo p ) — (ct/47t) ^jLo c /eD\\d c , and evaluating the above integral at zero temperature we 
obtain the zero-point-motion contribution 

( ^ )(0 ) = In ( i+V^ ^ (B7) 

in the case s > 1, which is physically relevant for the cuprates, and 



(f*J>(0) = _ ^ 2 arctan^^ (B8) 



in the case s < 1. For s — ► 1 both the above expressions (B7) and (|B8| ) reduce to 2uj c /tt£w p . For s ^> 1, from Eq. 
(B7), (5^)(0) ~ (8/ct) ln(^rWw c /(e-D||rfc)) in agreement with the expression used in the text, except the relevant 



energy scale, w c = 47re 2 /d c , for the anisotropic Coulomb potential appears instead of 27re 2 /£o, the scale appropriate 
when large Q.. contributes and we are closer to a 2D limit in the fluctuation integral. The quantum fluctuations are 
large for small s and decrease monotonically with increasing s. For s — > 00, (SO 2 ) — > 0, and the classical limit is 
recovered. 

At finite temperature and for a bare stiffness which is independent of temperature, the analytical result for the 
quadratic temperature dependence of (89 2 )(T) for arbitrary a, given in Eq.(|10|), has been derived in the text. The 



same result also follows from a low temperature analysis of Eq.(B6) above. 
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FIG. 1. Low temperature behavior of the renormalized superfluid stiffness for the dimensionless dissipation a = 150. The bare 
stiffness was chosen to be D ®d c = 80meV, independent of temperature, and the renormalized D^d c (T = 0) ~ 75meV, corresponding to 
(<5# 2 (T = 0)) ~ 0.1. The values obtained from the numerics are shown as squares while the solid line is the analytical T 2 form given in 
Eq. (|ic|). Linear temperature dependence sets in above a temperature ~ 20K as seen from the asymptote in the inset, which compares 
well with the estimate of the thermal crossover scale 3-D,, d c /a Kb 18K. 
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FIG. 2. The bare and renormalized 1/A 2 (T) plotted using dashed and solid line respectively, for a dimensionless dissipation <x = 150. 
The bare value of Ay and its slope were chosen such that the renormalized values, calculated using the SCHA equations (Q) and (^), 
Ay fs; 2100A and its slope dX^/dT fs; lO.OX/K, are in agreement with experiments Jli| ] on Bi2212. The renormalization due to quantum 
fluctuations is seen to be fs: 5%, much smaller than the ?s 50% renormalization obtained in the non-dissipative case IITJ] . 
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FIG. 3. The bare and renormalized stiffness plotted as a function of temperature, using dashed and solid lines respectively, to obtain 
an estimate of the transition temperature T c in Bi2212. We chose the bare bilayer stiffness d c D®(T) = d c D®(0) — oPT with parameters 
d c D®(0) ~ 80meV and a fs O.JmeV / K , relevant to bilayer Si2212. This leads to a renormalized bilayer stiffness of ss TomeV and a low 
temperature slope 0.7meV/ K for the bilayer stiffness, consistent with low temperature penetration depth experiments jl9| in _Bi2212. 
The renormalized stiffness is computed using the SCHA equations (^) and (^) with a dimensionless dissipation a = 150. The renormalized 
stiffness (solid line), within the SCHA, shows a jump near T ~ 90-R", but that is likely an artifact of the approximation. The temperature 
at which transverse excitations drive the transition to a non-superconducting state is estimated from the Nelson-Kosterlitz condition, as 
the point at which the dotted line intersects the renormalized stiffness curve (see text for details). This temperature, Tkt ~ 80A", is in 
reasonable agreement with experimental values [h9(| for T c ~ 90K in optimal Bi2212. 
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